Optically-controlled logic gates for two spin qubits in vertically-coupled quantum dots 
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We describe an interaction mechanism between electron spins in a vertically-stacked double quan- 
tum dot that can be used for controlled two-qubit operations. This interaction is mediated by 
excitons confined within, and delocalized over, the double dot. We show that gates equivalent to 
the VSWAP gate can be obtained in times much less than the exciton relaxation time and that the 
negative effects of hole-mixing and spontaneous emission do not seriously affect these results. 
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The spin of an electron confined in a quantum dot 
(QD) is one of the leading candidates for the realiza- 
tion of a practical qubit. Since the work of Loss and 
DiVincenzo there have been a number of proposals 
on how best to achieve the precise manipulations of such 
spins required for the operation of quantum logic. See, 
for example, @,&S@. 

Whilst interest in electrostatic gating remains strong, 
the use of lasers has several advantages in this role, most 
notably speed and control. Despite significant theoretical 
advances in this direction, there has, as yet, been no ex- 
perimental demonstration of optically-controlled gating 
between electron spins in QDs. 

In this paper, we describe an interaction mechanism 
to achieve just this. This qubit-qubit interaction is me- 
diated by intcrdot tunnelling of photo-excited carriers — 
an area which has been the subject of significant recent 
experimental advances 0, 0| ■ Our results are of explicit 
relevance to the current generation of vertically-stacked 
self- assembled InAs QDs, but are also easily adaptable 
to the other dots, including horizontally-coupled ones. 

The interaction we describe has its origin in the 
so-called optical- RKKY effect, in which two electron 
spins are coupled via their exchange interactions with 
optically-generated excitons in the semiconductor bulk 
[g]. The coupling effect of these bulk excitons between 
two electron spins in a double QD was examined in Ref. 
Q. Here we consider an interaction mediated, not by 
bulk excitons, but by a single exciton confined in the 
same double QD structure. We describe a situation in 
which the excitonic electron is able to tunnel between 
the dots and form delocalized 'molecular' states. It is 
the exchange interaction between this electron and the 
resident qubit electrons that leads to an optically con- 
trolled gating. This gate, although not one of the stan- 
dard quantum computation (QC) gates, can be used to 
form a controlled- operation when used twice in con- 
junction with single qubit rotations, and is, in this way, 
similar to the VSWAP gate. 

The main factor limiting the speed with which oper- 
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FIG. 1: (colour online) (a) Sketch of the two vertically- 
stacked quantum dots. Growth direction is the z direction 
and the two dots are positioned symmetrically at z = ±d. 
The dots have heights 2a and the width of the barrier between 
them is b = 2(d — a). The lower (A) and upper (B) dots have 
Darwin radii of La and Lb respectively, (b) Two electronic 
and two hole levels in each dot take part in gate operation. 
The lowest electron levels in each dot (labeled 1 and 2) are the 
two qubit levels and an electron permanently resides in each 
of them. Laser illumination is tuned such that it creates an 
exciton in the excited levels of dot A only (levels 3). (c) Due 
to a tunnel coupling T34 between the dots and a resonance 
condition met through the tuning of the gate voltage V g , the 
excitonic electron can tunnel back-and-forth. The exchange 
interaction between this electron and the qubits gives rise to 
the optically-induced interaction between the qubits that can 
be used to perform a quantum gate. 



ations can be performed in this set-up comes from the 
'kinetic exchange' between qubit spins which arises from 
the virtual tunneling of the qubit electrons between dots. 
This is an important consideration in vertically-stacked 
dots, and we show here that by choosing appropriate dot 
parameters, we can make this effect small and still obtain 
fast two-qubit operations. 

Our proposal therefore offers an accessible path to the 
demonstration of quantum logic in the solid state. Fur- 
thermore, it offers insight into the dynamics of interact- 
ing few-body systems in confined nanostructures — a 
topic of increasing experimental relevance. 

Figure [1] illustrates the basic principle at work here. 
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The two qubit electrons reside, one a piece, in the ground- 
state levels of a vertical double-dot structure. The dots 
have slightly different diameters and a static gate volt- 
age V g is applied in the growth direction. The distance 
between the dots is small enough that tunneling can oc- 
cur between the dots, but large enough to suppress the 
direct exchange interaction between the two qubit elec- 
trons. The tunneling of the qubit electrons into each 
others dot is strongly suppressed by the energy required 
to doubly charge one of the dots. 

A circularly-polarized laser pulse is used to generate 
an exciton into the lowest excited orbitals of dot A as 
shown in Fig.Qi>. The gate voltage is tuned such that the 
energy of this three-electron-and-one-hole configuration 
is resonant with that shown in Fig. [Ik, where the excited 
electron is in dot B. 

The resonant tunneling of the electron between these 
configurations creates a set of hybridised states in which 
the excitonic electron is delocalized over both dots. The 
qubit electrons experience a strong intradot exchange in- 
teraction with this delocalized electron and this medi- 
ates an effective interaction between the two qubits. By 
controlling the laser parameters, this interaction can be 
harnessed to perform quantum logic operations. 

We begin our detailed treatment by describing the 
single-particle states of the double dot. This includes a 
calculation of the size of the interdot tunneling elements. 
For reasonable dot separations we describe a situation 
where the qubits in the ground states are effectively iso- 
lated from one another, but where the exciton electron 
is free to tunnel. We then detail the spectrum of the 
four-body system resulting from the interplay of confine- 
ment, tunneling and the Coulomb interaction. From this 
follows the quantum logic properties of the system. Ana- 
lytically, we derive an effective gate operator for the sys- 
tem under adiabatic conditions. We also use numerics 
to explore the non-adiabatic capabilities. We show that 
gates equivalent to the VSWAP gate can be obtained 
in times much less than the exciton relaxation time and 
that the negative effects of hole-mixing and spontaneous 
emission do not seriously affect these results. 



SINGLE-PARTICLE STATES 



plane corresponding to dots A and B: 

M _ f 1/2 m e u> 2 A p 2 z<0, 
v pW-\ 1/2 m e uj 2 B p 2 z>0. 



(2) 



Here, p = \/ x 2 + y 2 is the radial coordinate, m e is effec- 



tive mass of the CB electron, and ftiux = 
the confinement energy of dot X = A, B 
nian of the electron is then 



H 2 / (m e L 2 x ) is 
The Hamilto- 
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V z {z) + V p {v) 



Eg. 



(3) 



where E g is the band-gap. 

We will consider the two dots to be separated such 
that a variational treatment in terms of wave functions 
localized in each of the individual dots is appropriate. 
Let us label with 1 and 3 the lowest two levels in dot A, 
and with 2 and 4 those in dot B. Our starting point, 
then, is the set of four single-particle wave functions 
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x(z + d)r) s {w A ;p,0), 
X(z - d)r] s (u) B ; p, 6), 
X(z + d)T] p (uj A ; p, 6), 
\(z - d)r) p (u B ;p, 6). 



(4) 



In the growth direction, xi z T d) are square-well eigen- 
functions centred at z = ±d. In the xy plane, the func- 
tions t/ s (lux; p,0) and rj p (ujx',p, 0) describe the ground 
and first-excited Fock-Darwin wave functions for confine- 
ment energy u>x- These wave functions are given explic- 
itly in Appendix |XJ This set of wave functions is not 
orthogonal since, although the overlaps S12 = (4>i\4>2) 
and S34 = (03 1 04) are small, they are finite. We thus 
use their orthogonalized counterparts ipi( r ) as outlined 
in the Appendix. 

The evaluation of the Hamiltonian of Eq. ([3]) in this 
orthogonal basis gives us our approximate single CB- 
electron Hamiltonian. Correct to first order in the over- 
laps 5*i2 ~ S34, we have 
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(5) 



The potential experienced by a conduction band (CB) 
electron in the double-dot structure may be approxi- 
mated as follows. In the growth direction z we have a 
double square-well potential with wells of width 2a and 
depth Vq centred at ±d: 



V z {z) 



\z±d\<a, 
Vq otherwise. 



(1) 



The width of the barrier between the dots (the parameter 
often quoted in experiments) is given by b = 2{d — a). 
Confinement in the x-y plane is assumed harmonic with 
different confinement energies above and below the z = 



The diagonal energies are given by 
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hio A ; E 2 
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E g + 2hui A ; £4 = E z + Eg + 2hu>B, (6) 



where E z is the energy resulting from the z confinement 
and is the same for both dots. 

The quantities T12 and T34 are ground- and excited- 
state tunneling amplitudes. They may be approximated 
as 



T12 ~ T34 ~ 4:qVqBC 
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FIG. 2: Logarithm of T34, the tunneling amplitude between 
excited states in different dots as a function of the interdot 
distance d. The solid black line shows the approximate result 
of Eq. and the black points show a numerical evaluations of 
full single-particle matrix element. The dashed line and grey 
boxes show the analogous quantities for heavy-hole tunnel- 
ing. The dot height is 2a = 2 nm and the radial confinement 
energies were fiw\ = 35 meV and hio% = 30 meV for all car- 
riers. The vertical confinement depth was Vq = 680 meV for 
electrons and Vq 1 = 100 meV for holes. We took the electron 
mass to be 0.04mo in InAs and 0.067mo in the surrounding 
GaAs matrix. Heavy holes were masses m/,, = 0.34mo and 
mux = 0.04mo for both InAs and GaAs. 



where /3, B, and C are functions of Vq determined by 
the continuity of the z-dircction wave function (see ap- 
pendix). The dependence of T34 on the interdot distance 
is illustrated in Fig. [2] for typical parameters. Current 
experiments have largely operated at small interdot spac- 
ings where the tunneling is large. For example, Ortncr 
and coworkers [|| measured a value of T12 ~ 13 meV for 
d w 3.5 nm, which is in good agreement with the results 
plotted in Fig. O For our purposes, we require a much 
smaller tunneling such that localized states are still well 
defined. We will consider separations in the range of 
d ~ 8 nm (b « 14 nm), where the tunneling amplitude 
lies in the range 10-100 /xeV. 

In Eq. © we have also taken into account the effects 
of the gate voltage applied in the z direction by incorpo- 
rating a first-order shift in the energy levels of dot B by 
an amount V g Q- 

In InAs dots, the valence-band states have a predom- 
inantly heavy-hole (HH) character. The single-particle 
Hamiltonian for such holes follows in exactly the same 
way as above except that the parameters differ. In par- 
ticular, the HH mass in the growth direction, m^z, is 
different to that in the xy plane rrihj_. The heavier mass 
of the HH in the z direction, as compared with the elec- 
tron, means that the HHs are more localized than the 
electrons. This in turn means that the HH tunneling 
amplitudes between the dots are approximately one or- 
der of magnitude smaller than those for the electrons, as 
illustrated in Fig. O 



II. LIGHT-MATTER INTERACTION 

The multi-particle basis states for our problem are cho- 
sen by filling the appropriate single-particle levels. Prior 
to the application of the laser, there are two electrons 
in the system and these reside in the ground-state lev- 
els of the dots. Interaction with the laser generates an 
additional electron-hole pair in the excited-levels of the 
double-dot system. 

We consider here illumination with a <x + circularly po- 
larized laser propagating in the growth direction. Con- 
servation of angular momentum means that this excites 
an exciton consisting of a spin-down electron and a spin- 
up heavy hole. Here, we have neglected hole-mixing and 
assumed that the in-plane magnetic field is zero — points 
which we return to later. The laser is tuned such that it 
creates the exciton in the excited levels of the dots and 
not in the ground-state levels where the qubit electrons 
reside. In principle, illumination can create an exciton 
with a hole in either dot 

The block structure of the Hamiltonian for the inter- 
acting system reads 



Hi 



total 



rph 



T& t 
J X4 





V fi^ 
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n< 3 ) 



G T, 



rpe t 
1 G 



G 



(8) 



D ) 



The two blocks in the bottom righthand corner, G and 
D, describe states consisting of the two resident elec- 
trons only — in G the electrons are separated one in 
each dot; in D, one of the dots is doubly occupied. These 
two blocks are coupled via the electron-tunneling terms 
of block Tg. The blocks and X^ describe states 
with two resident electrons plus the exciton, with the 
superscript referring to the dot-level in which the hole 
is located. These two blocks are coupled by the hole- 
tunneling terms of block T h . Block X D describes states 
where an exciton is present and one of the ground-state 
levels is doubly occupied. This block is coupled to the 
other exciton-containing states through the electron tun- 
neling blocks Tf 34 . Finally, the blocks fiW; { = 3,4,D 
contain the light-matter interaction terms, and provide 
the connexion between the ground and excited sectors. It 
should be noted that the effects of Coulomb interaction 
are incorporated into the diagonal blocks of iftotaii and 
that the exciton blocks, X« also contain single-electron 
tunneling terms. 

We can make several approximation that reduce the 
complexity of fftotai significantly. Firstly we assume that, 
whereas in principle the laser can excite an electron hole- 
pair in either of the dots, we assume that the two dots 
are sufficiently distinct in size that, in fact they are only 
generated in dot A. We will later chose the difference in 
confinement energies of the two dots to be w 5 meV, so 
this is certainly valid. We can also neglect the generation 
of the indirect excitons such as e^/i^, as these have very 
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small transition matrix elements. This means that the 
light-matter interaction is well approximated as 

H int = fi^e-^eJA + H - c -< ( 9 ) 

where Q(t) is the time-dependent Rabi frequency of the 
exciton transition. The block Q ( 4 \ corresponding to the 
generation of excitons in dot B, is therefore zero. Fur- 
thermore, hole-tunneling is very effectively suppressed in 
this system — not only are the hole-tunnelling matrix 
elements much smaller than their electronic counterparts 
(see previous section), but under the operating conditions 
of our device, to be described below, these transitions are 
far off-resonant. We therefore treat the block X^ 4 \ with 
a hole in dot B, as being effectively decoupled from the 
dynamics that we are interested in. 

Finally, transitions to states in which the dot lev- 
els are doubly occupied are suppressed by the enhanced 
Coulomb interaction of such states. That this is still the 
case even with different dot sizes and in the presence of 
applied field is jusitified in the next section. This means 
that the blocks D and X D are decoupled from the rest of 
the system over the relevant timescales. Taken together, 
these approximations mean that the system can be de- 
scribed with just the portion of the Hamiltonian -fftotai 
consisting of the blocks G, X^ and Since these 

blocks have significant internal structure themselves, the 
solution of this problem is still by no means trivial. 

III. GROUND-STATE SECTOR: ISOLATED 
QUBITS 

We continue our treatment by temporarily reinstating 
block D and consider all states in which there are only 
two electrons in the system. This is the "ground-state 
sector" , and is described by blocks G, D and Tg above. It 
is clear that only states with one electron in each dot rep- 
resent good two-qubit states and that double-occupancy 



of either dot, as in block D, is a source of error. It is 
therefore the purpose of this section to describe this sec- 
tor and to show the extent to which block D can be 
neglected and the qubits thought of as isolated. 

We can decompose the proper qubit sector as G = 
Diag{G_i, Go, G+i}, where the subblocks are labelled 
by the z-projection of the total electron spin. Block 
G-i corresponds to the energy of the two-qubit state 
ej^e^jO), G + i corresponds to e^eJ^lO), and Go consists 
of the two qubit states with the electrons having opposite 
spins: e^e^lO) and ej^e^jO). 

The exclusion principle determines that there are only 
two states in block D, namely e[,e|||0) and e^e^lO). 
Conservation of electron spin angular momentum means 
that only block Go couples to the block D. We can there- 
fore consider the effects of double occupation by concen- 
trating on the following Hamiltonian 

each block of which is a two-by-two-matrix. 

To assess the effects of double occupancy it is neces- 
sary to consider the Coulomb interaction between the 
confined carriers, in this case electrons. The electron- 
electron matrix elements between single-particle orbitals 
are given by 

Vm =Ec J d Tld ra |n-r 2 | ^ 

with Ec = e 2 /(47reo£i?). We define similar quantities 
Vijkl f° r electron-hole interactions and will only consider 
electron-hole direct interaction as e-h exchange effects are 
negligible. 

Including both Coulomb interactions and single- 
particle terms, the Hamiltonian of the ground-state sec- 
tor in question is 



V a 



H GD 



Jl2 

V -rf 2 



J 



12 



E l 



V a 



'12 



E l 



ceec 
221 



Vfffl - AE 







rpG T/eeee 



'12 



T/eeee 
^2222 



(12) 



AE + 2V g J 



where we have defined E G = Ef + E% + V^ff > = For the dot separations in which we are interested here 

E% — Ef, and J 12 = Vf^ii is the magnitude of ground- (d w 8 nm), this exchange energy is of the order 

state exchange interaction. This exchange term can be 10 -5 mcV and can thus be neglected, 

approximated as The tunneling amplitude Ti2 means that double- 
occupancy of the qubit levels can, in principle, occur, 

_ 4 and this would clearly be deleterious. However, if the 

J 12 fa cyx a ^ e -i/3d^ ^3) doubly occupied states are much higher in energy than 

L 2 A L 2 B {L^ + L^ 2 ) 3 ^ 2 the singly-occupied ones, double occupancy will be sup- 
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pressed . This requires that the energy differences satisfy 

Qi = Vini ~ ~AE-V g » rf 2 , 

Q 2 = V£H - V^t + AE + V g » r{ 2 . (14) 

With these inequalities satisfied, the effect of double oc- 
cupancy is reduced to imparting a small residual 'kinetic 
exchange' interaction on the fixed qubits. To second or- 
der in T12 this kinetic exchange energy may be approxi- 
mated as t 2 2 /Q, with Q^ 1 = Q^ 1 + Q 2 l ■ The important 
time-scale, 

Tq = hQ/T 2 l2l (15) 

follows accordingly such that, for times t <C Tq, the 
qubits are effectively isolated from one another. 

We will see later that device operation requires that 
the offset V g should be set as 

j-ie j-ie | Treeee Treeee , yreeee ireeee 

g — ^3 — ^4 "T ^1331 — v 2442 T ^2332 — v 144,1 

+ ^4334 G ~ ^3333 e + ^24 - ^13- (16) 

As the Coulomb elements are approximately symmetric 
between the dots, we can approximate Qi as 

Q_ y/eeee \/-eeee , Trehhe \/-ehhe , a 77 

1 ~ V U11 - V \22\ + ^3333 _ ^4334 + Aii , 

r^j T/eeee x/eeee rrehhe , j/ehhe a p (~\ *7\ 

W2 ~ V 2 222 ~ l/ 1221 - ^3333 + M334 _ I 1 ' j 

In order to obtain as long a time Tq as possible, we 
need to maximize Q, and this is achieved when Q\ = 
Q 2 . This in turn means that the optimal difference be- 
tween the confinement energies of the two dots should 



be AE ks -VfJ$. e + Vg^\ e , in which case Q = Qi/2 w 
l/2(Vf^f - VfHf). Typical parameters give AE w Q « 
5 meV and thus Tq = 1.3 ns for a tunneling amplitude 
of T12 = 0.05 meV, and Tq = 82 ps for t 12 = 0.2 meV. 
Providing then that the operation of our gate is much 
faster than Tq, we can neglect the effects of double oc- 
cupancy in the ground-state sector. A similar analysis 
shows that tunneling of ground-state electrons is sup- 
pressed to a similar degree also in the presence of the 
exciton. This means that, with the above proviso, we are 
justified in neglecting blocks D and X D from H to ta,\- 

The other important time-scale for this system is the 
relaxation rate of the exciton. Experiments have shown 
that this quantity is of the order of T le \ us 1 ns or greater 
[13, EH- This, therefore, is of the same order of magni- 
tude as Tq for t\ 2 — 0.05 meV. 



IV. EXCITON STATES 

We now consider states with the exciton, and since 
the hole is always localized in state 3, we consider only 
block X<- 3 \ We can decompose this block into subsec- 
tors based on the z-projection of the total electron spin: 

= mag{x^ 2 ,x!tf 2 ,X { *l /2 ,X { *l /2 y Since there 

are no spin-flip processes effective over the time scale of 
a typical gate operation, these sectors are independent 
and each excitonic sector couples exclusively to one of 
sub-blocks of G. This means that the total Hamiltonian 
for the system can be written 
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G+i 





^ o 

















y(3) 
^+3/2 



(18) 



Let us label each of these interacting sub- 
blocks with the electron-spin projection of the 
exciton state. We therefore have H to ta\ = 

Diag{ J fr(- 3 / 2 ) ! F(- 1 / 2 ),ir( 1 /2) )jH -(3/2)} Block #0/2) is 

clearly decoupled from the dynamics of the system and 
can be ignored. 

Let us consider the sector described by i/(~ 3 / 2 ', in 
which both qubits are spin-down and thus parallel to 
the electronic spin of the exciton. Only three states 



then make up this sector: the photo-excited state |I) = 
e il e 2i e 3i I ' the same configuration but with the elec- 
tron tunnelled |II) = e^e^e^/i^O), and the qubit 
ground-state | III) = e^e^jO). 



6 



The diagonal elements of the Hamiltonian for the first 
two states are 

h ujT ] = C u + 2V g -J 2 4, (19) 

where C\ is the sum of all single-particle energies and 
direct Coulomb interactions for state \T); 

d = E{ + El + El + El 

iT^eeee _i_ i^eeee _i_ Greece 
+ ^1221 v 1331 ' / 2332 

jrehhe Trehhe j/ehhe /OPA 

v 1331 ,/ 2332 v 3333 ' l zu 7 

and similarly for Cn. The energies J13 = V^ffl and J 2 4 = 
^2424 are the intradot exchange interaction strengths, 
and all other exchange interactions are negligible. 

In order that the excess electron can tunnel between 
the two configurations I and II, we set the voltage V g 
across the device such that we align = H^ u = E x . 
This requires 

V g = C\ — C\\ — J13 + J24 

rpe Trie 1 r/eeee \reeee < i/eeee T/eeee 

— H/ 3 Hi i -f- K 1331 ^2442 T ^2332 v 14il 

1 -irehhe j/ehhe , 7 7 /ni \ 

+ * / 4334 - K 3333 + ^24 ~ ^13, (/J-J 

which yields V g ~ 5 meV for our parameters. 

With this condition, the Hamiltonian for this sector 
reads 

/ E x -t 34 ne~ iuJt \ 
= -r 34 E x • (22) 
\ VLe iut E G J 

Moving to a basis of tunneling eigenstates for the exciton 
levels, 2" 1 / 2 (|I) ± |II)), we have 

( E x - t 34 Qe~ iut /V2 \ 
i? ( ^ 3/2) = E x +r 3i Cle- iut /V2 -(23) 

yf2e iw VV2 fie M 7V2 £ G / 

We tune the laser frequency such that, apart from a small 
detuning 5, it is on resonance with the transition from the 
ground state and the lowest tunnel coupled exciton state 
with energy E x - t 34 : 

uj + 5 = E x -T34- E G . (24) 

Finally, we move to a rotating frame 

H i™) H (m) = RH (m) R i + iRR ^ ( 2 5) 



with, in this case, R = Diag {e^+Eo)^ e i( - u+Ea ^, e lEot } 
to obtain 

/ S 0/V2\ 

H ( R 3/2) = 5 + 2r 34 0/V2 • (26) 

\ n*/V2 n*/V2 7 

We now consider the Hamiltonian, \ for the 

sector in which the qubit spins have opposite directions. 
This is the important sector as it is here that the ex- 
change between exciton and qubit electrons is manifest. 
Taking into account both electron tunneling and electron 
exchange, there are six exciton states in this sector. With 
intra-dot exchanges J13 and J24 diagonal, these are 

|I} = e^e^lO), 
|II) = 2-Va( e t T 4 i 4 i + e t ie t ie t T )^| 0)j 

[HI) = e^e^e^O), 

|IV) = 2- 1 /2^ e t ;e t Te t i + e t ue t ie t T )^| 0)) 

|V) = 2 I [e\^ 2 ^e\^ — ej^e^e^ ^.JO), 
|VI) = 2- 1 / 2 (e t ue t T ei i -e t u et i el T )^|0). (27) 

The first four states are degenerate under the intradot 
exchanges, and states 5 and 6 are split from the rest 
by an energy of 2J13 w 2J 2 4 ~ 6.5 mcV. This means 
that they are effectively decoupled and can be neglected 
henceforth. 

We can then proceed to the basis of tunneling cigen- 
basis defined by 

|I') = 1/V6|I) + 1/V3|II) + 1/V6|III) + 1/V3|IV), 
|II') = 1/V3|I)-1/V6|II)-1/V3|III) + 1/V6|IV), 
|III') = -l/V3|I) + l/V6|n)-l/\/3|III) + l/\/6|IV), 
|IV') = -l/V6|I)-l/V3|II) + l/V6|III) + l/\/3|IV). 

(28) 

In the basis of these four states plus the two ground states 
e^eJji and e^e^ the Hamiltonian in the rotating frame 
reads 
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fi*(t)/\/l2 -0*(t)/\/6 



n(t)/V6 
n(t)/Vs 

-0(i)/V3 
-fi(i)/\/6 







0(i)/V6 \ 
-Sl(t)/Vl2 
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(29) 



where the rotating frame is defined with frequency uj+Eq 
for all excited states and Eg for the two ground states. 

Analysis of the final sector, in which both qubit spins 
point up, proceeds exactly as for i?( -1 / 2 ) above, except 
with all electron spins flipped and only a single ground 
state e^e^O). The resulting Hamiltonian, H^/ 2 \ is the 
same as that of Eq. (f2"9"]) but with the fifth state omitted. 

Taking these results together, and ignoring for the mo- 
ment the light-matter coupling, we see that the excitonic 
sector consists of four energy levels, of which the inside 
and outside pairs are doubly and triply degenerate, re- 
spectively. In each of these states the excitonic electron is 
delocalized over the entire structure, and is thus capable 
of mediating an exchange interaction between the qubits. 
Note that the optical coupling between these states and 
the qubit ground states shows a variety of different cou- 
pling strengths, as determined by the forefactors of f2 in 
the preceding Hamiltonians. 



V. QUANTUM GATES 

Having derived these Hamiltonians, we now show how 
the interactions they entail may be used to perform two- 
qubit operations. We begin by considering the system 
in the adiabatic limit in which the exciton levels are but 
virtually populated. This approach yields an effective 
gate operator U c r acting only on the qubit space. Due 
to the complexity of the exciton structure, U c r does not 
have the form of one of the standard QC gates. However, 
we give the relationship between J7 e ff and the controlled- 
phase (CPHASE) gate, which makes clear the QC ca- 
pacity of U e g . We then consider non-adiabatic operation 
of the system and show how this can improve operation 
times. 

Second-order Rayleigh-Schrodinger perturbation the- 
ory of Raman processes can be used to derive an effective 
Hamiltonian for the qubit sector through the elimination 



of the exciton levels. In the basis {| J.J.),| |t), | ||),| TT)} 
we obtain the effective Hamiltonian 
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The validity of this Hamiltonian is conditioned on the 
standard adiabatic conditions of TS/H 3> 1, with T 
the operation duration, to avoid populating the trion 
levels, and |0(i)|/<5 -C 1 such that perturbation the- 
ory is valid. Under these conditions, the time-evolution 
operator of the qubit sector may be approximated as 
U cS (t) w er lAKl with A = - / dt\n(t)\ 2 /h8. We will use 
a Gaussian pulse shape, Q(t) = A exp (— t 2 /2T 2 ) . such 
that A = — A 2 y/irT/h5. Evaluating the matrix exponen- 
tial we find that C/ e rr can be written as 



e iq - cos 8 

-ie^+i+V 2 sin9 
\ 





_j e i(9-+?+)/2 sin 

e iq+ cos e o 

e iK < 



(32) 



8 



where Ki = —Axi, K3 = — Ax'3, sin 8 = sin <p cos 9 with 
4> = A/2y/ (x 2 — X3) 2 + Ay 1 and tan 9 = (x 2 —xz)/2y, and 
where q± = ±\ — A(x 2 + x^)/2 with tanx = sin 9 tan <f>. 

As is evident, U c g is not a standard QC gate. How- 
ever, we now make the connection between U e g and the 
familiar CPHASE gate 



C($) = Diag I 



,1,1,1 



(33) 



where $ is the controlled phase. Let us define the single- 
qubit operation Si (a) = Diag (e ia , l) acting on qubit i, 
and the generalized SWAP gate 
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(34) 



Setting 9 = tt/2 in W(Q) yields a SWAP gate, and set- 
ting 9 = 7r/4 gives a VSWAP entangling gate. In terms 
of these operators the physical gate U e g can be written 
as 



U, 



off 



<1- 



^)S 2 (<± 



K>3 



x C{<P/2)W{Q) 



xSi( 



9- 



^)S 2 ^ 



K 3 



). (35) 



where 
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(36) 



The gate U g therefore contains a product of two entan- 
gling operations, C($/2) and W(8). 

We can isolate the CPHASE gate from the product 
CW by observing that the key unitary transformation of 
W by the single-qubit operation S' 2 {ir) = Diag (l, e~ l7r ) 

gives S' 2 (tt)W S' 2 \tt) = . Since S 2 commutes with C, 
the concatenation of CW and its transform yields, 

C(<S>/2)W(Q)S 2 {Tr)C(<S>/2)W(Q)S' 2 (ir) = C($). (37) 

Note that the transformation is independent of 8 and 
can be carried out without the knowledge of its value. 
From this it follows that the composite operation of two 
applications of U e g interspersed with appropriate single- 
qubit operations is equal to a CPHASE gate with the 
controlled phase $ as defined in Eq. (|36l) . 

The entangling capability of U e g can be assessed by 
calculating the concurrence C of the states formed by 
acting with U c g on a separable state of qubits 1 and 2, 
and averaging over all such inputs. This we find to be 



C = l/8|15-8cos($/2)cos(28)-7cos(48)| 1/2 , which 
depends only on the angles $ and 8. The average concur- 
rence is thus bounded by \ |sin($/4)| < C < y/ll/32 w 
0.586, where the lower bound comes from the 8 = case. 
In the example that we consider later (operation 1 of Ta- 
ble |TJ), the off-diagonal angle is 8 = 0.42 and the average 
concurrence is C = 0.50. 

We will quantify the usefulness of gate U g for QC pur- 
poses solely through the CPHASE angle <&, even though 
the gate U really requires two angles $ and 8 to spec- 
ify its entangling action completely. This one-parameter 
characterisation is appropriate because, from the forego- 
ing, we know that knowledge of $ sets a lower bound 
on the average entangemcnt generated by the gate. Fur- 
thermore, via the construction of Eq. (|37p . we know ex- 
plicitly how to form a full CPHASE gate independent of 
the value of 8, and indeed independent of the knowledge 
of its value. Clearly for U c g to be an efficient gate we 
require 4> to be significant. A value of $ = tt gives the 
composite operation as a controlled-Z operation, and in 
which case the constituent U e g operator may be said to 
be 'as efficient' as a VSWAP gate, since two applications 
of vSWAP are also required to form the controllcd-Z 
operation [![. Experimental demonstration of U e g with 
$ = tt would, in this sense, be equivalent to the demon- 
stration of a %/SWAP gate. 

The question then becomes whether we can obtain a 
gate U c g with a value of $ = ir within the time con- 
straints set by Tq and Tr, c i. The answer is that if we 
insist that the evolution of the system be strictly adia- 
batic, then for a typical tunneling of T34 = 0.05 meV, 
such a gate takes several hundred picoseconds. At this 
coupling, we have Tq = 1.3 ns, and so the adiabatic pulse 
time T is a significant fraction of Tq, which is undesir- 
able. However, if we allow a degree of non-adiabaticity 
in the time evolution, we can significantly reduce this 
operation time, as we now discuss. 

If we populate the exciton levels then we require that 
the population is returned to the qubit sector at the end 
of the operation. We quantify the degree to which this 
happens with the average norm on the qubit space. Let 
us write as [/total the time-evolution operator of the total 
system. The corresponding operator for the qubit space 
is U = Pf/total-P, where P is the projector onto the 4x4 
space. We can then define the average norm on the qubit 
space as No = j Yli i=i Wij | 2 - A perfect return has N = 
1 and in this case U necessarily has the same form as 
Eq. l[3"2"j). given the Hamiltonian of Eq. ([15]). 

We then numerically search the space of operations 
possible within the system for valid gate operations. For 
a given tunneling strength T34, pulse duration T, am- 
plitude A, and detuning 8, we numerically integrate the 
time-dependent Schrodinger equation for the Hamilto- 
nian of Eq. (fT8|) to obtain the systems evolution from the 
four initial qubit states J, J,, |t> Tl an( l TT- From these re- 
sults we can construct U, the evolution operator on the 
qubit space. We consider that the time-evolution can be 
construed as a valid gate if Nq > 0.99. 
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T34 


T 


T/T Q 


5 


A 


$ 


No 




1 


0.05 


84 


0.064 


0.0729 


0.10983 


3.1415 


0.991 


0.970 


2 


0.05 


125 


0.095 


0.084 


0.078025 


3.1415 


0.99993 


0.984 



TABLE I: Parameters and results for two example operations. 
T34 is the tunneling amplitude between dots, and A, T, and 5 
are the amplitude, duration and detuning of the laser pulse. 
The units of T34, 5 and A are meV, and T is measured in ps. 
For both operations $ ~ tt, as require. Without spontaneous 
emission the fidelity is equal to the norm on the qubit space 
T = No and this can be made very close to unity, as in the 
second operation. In the presence of spontaneous emission, 
the fidelity jT rol is slightly reduced. 



FIG. 3: (colour online) The absolute value of the optimal 
CPHASE angle $ pt as a function of T/Tq, the pulse dura- 
tion scaled by the qubit isolation time. Results are shown 
for three different tunnel strengths T34 = 0.05 meV (blue 
squares), T34 =0.1 meV(green triangles), and T34 = 0.2 meV 
(red diamonds). As is evident, a phase of $ = tt can be ob- 
tained at T 34 = 0.05 meV with T/Tq > 0.064. In real time 
units this equates to T « 84 ps. This value of $ can also be 
achieved with higher tunneling amplitudes, but this requires 
longer pulses. 



VI. HOLE-MIXING AND SPONTANEOUS 
EMISSION 

Other than the direct and kinetic exchanges between 
the qubits discussed above, there are two further dele- 
terious processes that we now discuss: hole-mixing and 
spontaneous emission. As we know from the Luttinger- 
Kohn model of semiconductor bands [HI, EH, hole-mixing 
means that, rather than the 'bare' heavy-hole states 

I 5 
by 



± I ) , the actual hole states are better approximated 



For a fixed pulse duration T, we obtain a range of valid 
gate operations with various different values of A and S. 
From this set, we are interested in the gate with the most 
significant gating action which, in terms of the CPHASE 
angle described above, is the gate for which |$| is closest 
to 7r. We denote the CPHASE angle of the gate which 
fullfills the criteria at a given T as $ pt- 

In Fig. [3]we plot the optimal angle $ opt as a function of 
pulse duration T for several different tunnel amplitudes 
T34. Let us consider the r 34 = 0.05 meV case. The angle 
<& pt reaches a value of tt at a value of T/Tq = 0.064 
and then saturates. At this coupling, Tq = 1.3 ns and so 
the pulse duration is T — 84 ps, which is short compared 
with both Tq and T TQ \. For the higher couplings shown 
in Fig. [3l wc obtain $ op t = tt only for ratios T/Tq > 0.1. 

In Table HI we give parameters and results for two ex- 
ample T34 = 0.05 meV operations: the first corresponds 
to the minimum pulse time T required to obtain $ = tt, 
and the second is a longer duration pulse with higher 
norm N Q . In the course of the evolution with the first 
set of parameters, the exciton states are occupied with a 
peak population of the order of 40%. For initial condi- 
tions of three of the four basis states virtually no popu- 
lation is left in the trion levels at the conclusion of this 
operation. However, evolution from J,f does leave a small 
remnant and this reduces the qubit norm to its value of 
Nq = 0.991. For longer pulse durations, as in the second 
operation of Tab. [H this remnant can be eliminated. 

Permitting the system to evolve non-adiabatically then 
improves the performance of the gate by reducing the 
time required to perform a useful rotation. 



cos ( 



3+3 
2 ' 2 



sin tie 



, (38) 



where | | ,=F \ 



are light-hole states, and 6 m , (j) m are 
mixing angles. With this mixing, the light-matter inter- 
action Hamiltonian becomes 



VL{t)e- wt {cos0 r 



sin 



0™4i h zH-} +H.c.,(39) 



where the factor of \/l/3 in the second term comes 
from the different weights of in-plane components of the 
valence-band wave functions [HI, [H| • 

This interaction requires that we consider an extra 
excited-state sector with \H~) instead of \H + ). Tak- 
ing both sectors into account, we can readily derive an 
effective Hamiltonian in the presence of hole mixing. We 
find 



Wmix = cos 2 e m H cf f + ^ sin 2 6 m H- , 
where 7i G ff is as in Eq. ([30|) and 



(40) 



H- = — - 
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(41) 



with terms as in Eq. (|3 1 [) . Note that the hole-mixing does 
not introduce any new matrix elements interactions be- 
tween qubit states that were not already present in the 
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non-hole-mixing Hamiltonian. From this and the sym- 
metry of ?i m i x it immediately follows that the CPHASE 
angle is $ m i X = {cos 2 6 m + | sin 2 6> m } where <3? is the 
original angle of Eq. ([55]) . Thus, the effect of hole-mixing 
is to slightly decrease the CPHASE angle <!>, and this can 
simply be compensated for by a proportionate increase 
in Rabi frequency f2. 

We now consider spontaneous emission. Given our 
laser excitation and the magnitudes of the indirect oscil- 
lator strengths, the dominant spontaneous emission path 
is the direct recombination of exciton pair e^h^. As 
noted previously, the time-scale of this process is of the 
order of T le \ = 1 ns. We can assess the effects through 
the numerical integration of the master equation for the 
system in the Lindblad form. 

We judge the quality of the operation in the presence 
of spontaneous emission through the fidelity [3 HI, de- 
fined as 



T= <*in|C/VoutC/|*in), (42) 

where the overline represents an average over all input 
states l^in), Pout is the output density matrix given |\&i n ), 
and U is the ideal gate operation. We define U by tak- 
ing the gate operator without spontaneous emission and 
renormalising it such that Nq = 1. Correspondingly, the 
fidelity of the operation U without spontaneous emission 
is simply given by T = No- 

The results of these calculations can be appreciated by 
again studying the two operations listed in Table fl] In 
this table we give the fidelities in the presence of sponta- 
neous emission. For the first, shorter pulse, spontaneous 
emission reduces the fidelity from 0.99 to 0.97. For the 
longer pulse, fidelity is reduced from very close to unity 
to 0.984. This final value is typical for operations with 
spontaneous emission, with the final fidelity lying in the 
98—99% range. 

This reduction in fidelity arises largely from population 
being left in the trion levels, as attested by the fact that 
the norm on the qubit space is affected to a similar extent. 
This residual population occurs because the spontaneous 
emission removes the direct exciton e^h^ but leaves the 
indirect exciton e^h^ untouched. Since, after the laser 
pulse has passed, this indirect state has no route back to 
the ground state, it is effectively trapped. This popula- 
tion eventually returns to the ground state, through ci- 
ther interdot tunneling plus direct spontaneous emission 
or, to a lesser extent, indirect emission, but this occurs 
on a time scale greater than that considered here. This 
trapping effect, although ever present, is small enough 
that operations of sufficient quality can still be obtained. 



VII. SINGLE-QUBIT ROTATIONS 

An important feature of the current set-up is that it is 
compatible with single-qubit operations on the electron 
spins. As described in Ref. [l4| (see also Ref. [Hj]), single 



qubit operations can be performed through the excitation 
of a exciton in the ground-state levels of a single dot 
(e.g. the levels 1 in dot A). The addressability of the 
individual dots is assured here by their different Darwin 
radii. To obtain arbitrary rotations, a static magnetic 
field is required in the x-direction (perpendicular to the 
growth direction) , and it is important that the two-qubit 
operations function correctly in the presence of this field. 

In InAs dots, both the electron and hole g-factors are 
finite [It], EH , and this leads to a splitting of the single- 
dot trion levels and the linear polarization of the optical 
transitions [3, [III ■ Assuming that both g- factors have 
the same sign, we can perform our two-qubit operations 
with a H-polarized laser (as defined in Refs. [3, E3) 
rather than a er + -circularly polarized one. We then de- 
rive an effective Hamiltonian for this situation and find 

Hb = \{H cS +H-{5 + Y> b )}- (43) 

where the second term is the same as Eq. (|41[) but with 
the detuning 6 replaced by 6 + throughout. Here 
S g = (gl + g^)fj, B B x , with g% and g x the electron and 
hole in-plane ^-factors, B x the in-plane field, and the 
Bohr magneton. As with the hole-mixing, we see that the 
magnetic field does not add any new matrix elements be- 
tween qubit states, and thus the operation has the same 
character as before. Now, however, the CPHASE angle is 
given by <f> B = 1/2 {$ + $(<5 + £b)}, where the second 
term is the same as Eq. (|36|) but with 5 + Y.b replacing S 
throughout. The contribution of this term is expected to 
be small since, with ^-factors \g*\ = 0.46 and \g%\ = 0.29 
[HI, El, at B = 8 T, we have |£ B | w 0.4 meV which is 
larger than typical values of 5. In any case, these changes 
to $ can easily be compensated for by a trivial change in 
amplitude of the laser. 

VIII. CONCLUSIONS 

We have described a mechanism for obtaining optical- 
controlled quantum gates in quantum dots. This may be 
thought of as a confined ORKKY interaction in which 
the exchange between the two qubit spins is mediated 
by a set of exciton states delocalizcd over the double dot 
structure. 

For this gating to function it is imperative that the two 
dots be spaced an appropriate distance apart such that 
the tunnel coupling between them is small enough that 
the kinetic exchange can be neglected over the course of a 
gate operation, but large enough that laser pulses short 
enough compared with exciton relaxation time can be 
used. Our calculations show that a distance of d « 8 nm 
would be ideal. In this case, significant gates (equivalent 
to ySWAP) can be obtained at high fidelities with pulse 
durations of the order of lOOps. 

The speed of this operation is essentially limited by the 
occurrence of the kinetic exchange. Without this effect, 
the dots could be closer together, the tunneling matrix 
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elements larger, and hence the strength of the gating in- 
teraction much stronger. However, this exchange is un- 
avoidable in vertically-stacked QDs, since such dots only 
grow in a stack when the separation between the two 
dots is small fl9| . This problem, therefore, is not unique 
to our gating mechanism, but rather a feature of cou- 
pled QDs, and any proposal seeking to show controlled 
quantum operations in such a structure must take it into 
account. 

In this respect, horizontally-coupled quantum dots 
have an advantage since, as such dots are coupled in 
the xy plane, the overlap of the exited states is a fac- 
tor of \[2 greater than that of the ground states. This 
means that the strength of the kinetic exchange interac- 
tion is reduced by a factor of one-half relative to exciton- 
tunncling. This in turn means that pulse-durations half 
that described here can be used to obtain the same re- 
sults. 

We have also demonstrated that our two-qubit opera- 
tions function equally as well with an in-plane magnetic 
field. This is important because it shows that our set-up 
is compatible the single-qubit rotations of Ref. [l4|. In 
combination, this means that we have demonstrated that 
the laser excitation of excitonic states provides a realis- 
tic, experimentally-accessible protocol for the fast perfor- 
mance of a universal set of quantum operations on two 
electron-spin-qubits in a coupled double quantum dot. 

This work was supported by ARO/NSA-LPS. 



APPENDIX A: SINGLE PARTICLE WAVE 
FUNCTIONS 



The ground-state energy E z may be found through solu- 
tion of the equation 



atan(aa) = f3. 



(A4) 



In the xy plane we have the ground-state and first 
excited-state Fock-Darwin orbitals 



r) s (u>;p,9) = y-e"^ p , (A5) 

V P ±(^P,9) = ^pe-^e^ 6 , (A6) 

with A = muj/h. Since the two p-orbitals are degener- 
ate, and angular momentum is conserved, we only need 
consider a single orbital and choose r\ v = rj p+ for con- 
creteness. 

The four complete wave functions 4>i(r) are given in 
Eq. ((4]), but these are nonorthogonal. As our variational 
basis we therefore use 



V'i(r) = 7i _1 {^l (r) -gi<fa(r)} , 

■02(r) = 7i _1 {0 2 (r) - 5i0i(r)} , 

^3(r) = 73" 1 {03(r) - g 3 <Mr)} , 

lfr 4 (r) = is" 1 IM?) ~ 93fo(r)} , (A7) 



Our single-particle variational wave functions are sep- 
arable into z and xy components. In the z direction, the 
wave functions are those of ground state of a square well 
of depth Vq and width 2a: 



Cexp(f3z) z < —a, 
X(z) = { Bcos(az) \z\ < a, 
C exp(— (iz) z > a. 



where 



a 



h 2 



P = 



2m e {V -E z ) 



H 2 



(Al) 



(A2) 
(A3) 



which are all orthogonal. Here, the overlaps Sij are de- 
fined Sij = J d 3 r <f>i(r)<f>j (r) , and the normalization coef- 
ficients are given by 



<n 

7i 

73 



l-Jl-S 2 2 /S 12 , 



l-2 gi S 12 + g 2 , 
l-2g 3 S 34 + gl 



(A8) 
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